High symmetry versus optical isotropy of a negative index metamaterial 
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Optically isotropic metamaterials (MMs) are required for the implementation of subwavelength 
imaging systems. At first glance one would expect that their design should be based on unit cells 
exhibiting a cubic symmetry being the highest crystal symmetry It is anticipated that this is a 
sufficient condition since it is usually assumed that light does not resolve the spatial details of MM 
but experiences the properties of an effective medium, which is then optically isotropic. In this 
work we challenge this assumption by analyzing the isofrequency surfaces of the dispersion relation 
of the split-cube in carcass (SCiC) negative index MM. We show that this MM is basically optically 
isotropic, but not in the spectral domain where it exhibits negative refraction. The primary goal of 
this contribution is to introduce a tool that allows to probe a MM against optical isotropy. 



PACS numbers: 78.20.Bh, 78.20.Ci, 41.20.Jb 

I. INTRODUCTION 

Driven by the desire and the opportunity to have opti- 
cal materials with tailored properties at hand, the field of 
metamaterials (MMs) attracts a steady increasing share 
of research interest. To loosely define the field, one may 
understand MMs as artificial structures made of sub- 
wavelength unit cells , called metaatoms. They predom- 
inantly affect the light propagation by a careful choice of 
their geometry and their arrangement; and not by the in- 
trinsic properties of the materials they are made of. MMs 
allow to tailor the flow of light well beyond what would be 
possible with naturally occurring materials by mimicking 
unprecedented optical properties. A large variety of new 
optical phenomena were predicted and experimentally 
proven on the base of MMs, where negative refraction 
attracted potentially the largest amount of interestP^. 
Metamaterials providing a negative effective refractive 
index are the essential ingredient to fabricate a perfect 
lens with an optical resolution limit well below Abbe's 
prediction 7 . However, and this is the major requirement 
currently not met by most of present MMs, this effective 
refractive index can only be meaningfully introduced if 
the dispersion relation is isotropic, i.e. for lossless media 
the isofrequency surface of the dispersion relation in 3D 
k-space must be spherical. Moreover, this isotropy has 
to hold also for evanescent waves, since they carry the 
subwavelength information for the desired super resolu- 
tion. To overcome this obstacle various approaches for 
obtaining isotropic MMs have been put forwarcPtSJ 

A systematic approach to design an isotropic magnet- 
ically active MM was proposed by Baena et al. 15. The 
first step is to choose a highly symmetric MM. Such MM 
is constructed based on a cubic unit cell arranged at a 
cubic lattice. It is evident that MMs with this high sym- 
metry behave isotropically in the quasi-static limit. Fre- 
quently the optical properties of such MMs have only 
been investigated for normal incidence (zero transverse 



wave vector) and then extrapolated towards finite trans- 
verse wave vectors, assuming isotropic medium response 
due to the symmetry of the structure. 

However, and this is crucial, even if the structure is 
operated in the subwavelength domain, the required op- 
tical isotropy is not straightforward, since the typical res- 
onance wavelengths, evoking magnetic effects, are com- 
parable to the structure size. There are only a few at- 
tempts where the MMs are probed by obliquely incident 
fields and for different polarizations, mostly restricted 
to transmission and reflection measurement J 16 " 2Q l Since 
even for an isotropic medium the reflection and trans- 
mission coefficients evidently depend on the angle of in- 
cidence, it is hardly possible to draw conclusion with re- 
spect to the isotropic behavior of the respective MM. To 
ultimately verify whether a highly symmetric metamate- 
rial behaves optically isotropic, it is necessary to calcu- 
late the dispersion relation uj = uo (k x ,k y ,k z ) where its 
iso-frequency surface k z = k z (fc x , k yi uo = const.) governs 
the diffraction and refraction properties of the MM 21 . 
Here, cj, fc x>y , and k z are the frequency, the transverse 
and the longitudinal wave vector components, respec- 
tively, where the latter is frequently termed propagation 
constant. The aim of this work consists just in evaluating 
this iso-frequency surface. 

For computing iso-frequency surfaces of a three- 
dimensional MM we take the complex permittivity of 
the material into account. This is crucial because most 
the band structure solvers assume lossless media which 
is certainly incorrect for metamaterials. We use instead 
a plane wave expansion technique that solves Maxwell's 
equations in the frequency domain for the periodic struc- 
ture. The technique solves the respective eigenvalue 
problem for the generally complex propagation constant 
of the Bloch modes 22 . Computing the propagation con- 
stant as a function of the transverse wave vector for a 
fixed frequency provides the iso-frequency surfaces. By 
using this method we specifically show for the split-cube 
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in carcass (SCiC)-MM with cubic symmetry that its op- 
tical response is isotropic only at low frequencies. By 
contrast this isotropy disappears in the spectral domain 
where the propagation constant, and hence the effective 
refractive index is negative. 

After having introduced the system under considera- 
tion in Sec. II, we show in Sec. Ill that the proposed 
MM fulfills all requirements usually imposed for its ho- 
mogenization. In Sec. IV we discuss the iso-frequency 
surface for the negative refraction regime and introduce 
a simple measure to quantify the optical isotropy across 
the entire sub-wavelength frequency domain. 



The SCiC was chosen because it has the highest pos- 
sible symmetry for a periodically arranged MM. In par- 
ticular, it is mirror-symmetric with respect to three or- 
thogonal axes. This excludes any effects resulting from 
first order spatial dispersion like chirality^ 15 * 20 * 31 *. Due 
to the cubic symmetry the three main propagation di- 
rections are equivalent. Hence, the optical response of 
the SCiC is supposed to be described by scalar effective 
material parameters in the quasistatic limit, i.e. the op- 
tical response should be isotropic. From these symmetry 
considerations one usually draws the conclusion that the 
SCiC might be an ideal candidate for an isotropic nega- 
tive index metamaterial. 



II. SPLIT-CUBE IN CARCASS DESIGN 

For our investigations we use the split-cube in carcass 
structure, which is a simplified version of the split-cube 
in cage structure 2 ^ containing less number of fine de- 
tails. Using the nested structures approach the SCiC 
is designed such, that one element of the unit cell that 
exhibits magnetic properties is inserted into another el- 
ement which shows dielectric resonances 2 ^. The SCiC 
unit cell consists of two silver parts embedded in silica 
(e = 2.25). Silver is regarded as Drude metal with a 
plasma frequency of cj p = 1.37 • 10 16 rad/s and collision 
frequency uj c = 8.5 • 10 13 1/P^. These parameters were 
adjusted to emulate the real experimental situation in 
simulations. 




FIG. 1: (color online) Schematic of the Split Cube in Carcass 
unit cell. The Carcass (a) serves to provide a plasma- like 
optical response resulting in a negative effective permittivity. 
The Split Cube cut (cut is shown for clarity) (b) provides an 
artificial magnetic response. 

The outer part, the Carcass, is a three-dimensional 
wire medium 26 ^ 27 ! and provides a dielectric resonance 
('negative permittivity' e). The inner part, the Split 
Cube, is a hollow cube with slits in the middle of the 
facets. It is the logical 3D extension of the symmetric 
split ring resonator^^S CO ncept and provides a mag- 
netic resonance ('negative permeability' ja). The details 
regarding the structure's sizes are indicated in Fig. [I] 



III. HOMOGENEITY 

To describe a MM as an effectively homogeneous 
medium, it is necessary to assure at first that light propa- 
gation within the MM is exclusively governed by a single 
Bloch mode. To be sure that this requirement is met, 
we also have to prove that the incident light field couples 
only to this particular mode. This is done by compar- 
ing the effective plane wave propagation constant, as re- 
trieved from reflection and transmission data of a finite 
structure, with that obtained from the dispersion relation 
of the lowest order Bloch mode 3 - 2 ^. 

We will determine the effective propagation const ant in 
the SCiC by using the S-parameter retrieval met ho d 3 - 3 -* 3 -^. 
The complex reflection and transmission coefficients are 
numerically calculated by the Fourier Modal MethocP^ 
where 31 x 31 Fourier orders were retained to achieve 
convergent results. By inversion of these scattering data 
the effective propagation constant k = k z and the effec- 
tive impedance Z can be determined. To certify that 
such an assignment of effective parameters to structures 
composed of only a small number of functional layers is 
meaningful and that the homogenization is valid, it is 
necessary to investigate the convergence of the parame- 
ters towards their bulk values 36 , therefore we determine 
the propagation constant of the Bloch eigenmodes of the 
infinite structure as well. 

This is done by calculating the T-matrix for a single 
period and solving the eigenvalue problem 

t (Jfj = exp(ifcA) (Jfj 

to obtain the propagation constant k = ^Bioch °f the 
Bloch modes, where A is the period in the main propa- 
gation direction. 

The results for the real and imaginary part of both 
propagation constants (finite homogeneous and infinite 
periodically structured media) and the formally intro- 
duced effective refractive indices n = ck/uj are shown in 
Fig. H 

For the periodic bulk material we get in general the 
propagation constants for an infinite number of Bloch 
modes (which remains nevertheless finite because of the 
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FIG. 2: (color online) Real (a) and imaginary (b) parts of 
the propagation constants and effective refractive indices (c) 
and (d) for the effectively homogeneous finite and the infinite 
periodic structures at normal incidence. The results for the 
finite structure (solid lines) are given for an increasing number 
of layers (1 layer - blue, 2 layers - green, 5 layers - red). The 
dotted black line corresponds to the values obtained from the 
Bloch mode with the smallest losses. 

numerical truncation of the number of plane waves re- 
tained in the plane wave expansion) but only the zeroth 
order mode, i.e. the one with the smallest imaginary 
part is shown. Clearly the values obtained from the fi- 
nite effectively homogeneous structure are converging as- 
tonishingly fast towards the values for the periodic bulk 
medium. Thus, two important conclusions can be drawn, 
firstly the homogenization of the periodic MM is feasible 
and secondly the optical response of a SCiC-MM consist- 
ing of only a few layers equals that of a bulk medium. 

The design purpose of the SCiC is a negative effective 
propagation constant, and thus an negative effective in- 
dex of refraction which is clearly achieved for frequencies 
around 170 THz (Figure of Merit -3?(n)/9?(n) = 0.35). 
Note that the effective permeability is dispersive but pos- 
itive in the investigated frequency range. Hence, the 
SCiC is a single negative MM resulting in the fairly small 
Figure of Merit. For frequencies less than 150 THz the 
SCiC is rather a strong absorber due to the large metal 
fraction. 

Since for the scattering problem at the finite system 
the coupling to different Bloch modes was rigorously con- 
sidered, we can conclude that the propagation of light 
through the structure as well as the coupling process is 
almost entirely dictated by the fundamental Bloch mode 
only, again underlining the validity of the homogeniza- 
tion procedure. This is clear as the propagation constant 
values for the finite and the infinite structures coincide. 
Although not shown here for the sake of brevity the prop- 
agation constants for oblique incidence are also rapidly 
converging towards the bulk values. Also the effective 
impedance is converging for normal incidence as well as 



for oblique incidence. 

Hence, we can fully rely on describing light propagation 
in terms of the dispersion relation assuming that light will 
couple to this fundamental mode only. Note that from 
the results above we can also conclude that the SCiC ful- 
fills all requirements for the homogenization of the MM, 
namely the structure is subwavelength (compared to the 
wavelength of the environment) and light propagation in- 
side the structure is determined by a single Bloch mode 
to which an external light field predominantly couples. 



IV. DISPERSION RELATION AND OPTICAL 
ISOTROPY 

The main property we are interested in is the optical 
isotropy of the SCiC, in particular in the spectral region 
around 170 THz where the propagation constant is neg- 
ative. To judge this we will proceed in calculating the 
iso-frequency surface. 

Due to the symmetry of the structure the complete iso- 
frequency surface can be constructed by only calculat- 
ing the dispersion relation in the irreducible Brillouin 
zone. The real and imaginary part of the propagation 
constant for a fixed frequency of 170 THz are shown in 
Fig. |3ja) and (b). For comparison we show in Fig. [3jc) 
and (d) the analytically calculated iso-frequency surface 
of an isotropic medium, which has the same refractive 
index as SCiC at normal incidence. 




k /k 
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FIG. 3: (color online) Real (a) and imaginary (b) part of the 
propagation constant of the lowest order Bloch mode in the 
first Brillouin zone. Real (c) and imaginary (d) part of the 
propagation constant of an isotropic medium with the same 
refractive index as the SCiC at normal incidence. 

At first we notice that both the real and the imagi- 
nary part of the propagation constant of the MM are not 
rotationally symmetric and monotonically changing with 
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increasing absolute value of the tangential component k t 
as for the isotropic medium. In particular, the depen- 
dency on the transverse wave vector is tremendously dif- 
ferent when compared to the isotropic medium. While 
the real part is monotonically increasing with k t in the 
TX-direction it is non-monotonously decreasing in the 
TM-direction. Also the imaginary part is strongly in- 
creasing in the TM-direction and only slowly varying in 
the TX-direction. Hence, any formally introduced effec- 
tive refractive index would explicitly depend on the trans- 
verse wave vector = (k xi k y ). Therefore, the intro- 
duction of a global effective refractive index is pointless 
since no additional information is obtained. Neverthe- 
less for paraxial wave propagation the introduction of a 
local effective refractive index is feasible. Near the T- 
point the isofrequency surface is approximately spherical 
where the validity of this approximation strongly depends 
on the frequency and the wavelength to cell size ratio as 
discussed later in detail. 

It should be mentioned that the choice of real valued 
k x and k y is arbitrary to a certain extent. In lossy me- 
dia like metamaterials, in general, the complex nature of 
the dispersion relation cannot be neglected. Hence the 
dispersion relation could also be calculated for complex 
valued k t . On the other hand it is pointless to provide 
these values as they are not accessible in any experiment. 
The tangential wave vector components are continuous at 
boundaries and only real valued plane wave solutions can 
exist in free space, therefore this choice reflects experi- 
mental constraints. 

To more quantitatively evaluate the optical isotropy of 
the MM, we monitor in the following the relative devia- 
tion between the numerically obtained propagation con- 
stant and the ideal spherical isofrequency surface. To 
quantify this deviation we assume that the effective re- 
fractive index at normal incidence is valid also for any 
other propagation direction. From symmetry consider- 
ations and the exemplary results in Fig. [3] we conclude 
that it is sufficient to investigate the dependency in the 
high symmetry TX and TM direction assuming that for 
a fixed value of |fc t | these points are extremes of the iso- 
frequency surface. In Fig. [4] the real part of the propaga- 
tion constant is shown as a function of the frequency and 
of the tangential wave vector component fct/&o> where 
fco is the free space propagation constant. The edge of 
the first Brillouin zone for the largest possible frequency 
(200THz) would be at k t /ko ~ 3 because of 

h _ 7T A(200Thz) 2 • 10 6 
k ) ~ak ) ~ 2a 200 -lO 12 ^ 

The most important frequency region here is the black 
domain of negative refraction where the Split Cube pro- 
vides the artificial magnetic response. Also the iso-error 
lines for the relative error 

^ _ \k z \ — | k z | ideal 
I (ideal 

with | fc z | ideal = ly^WWt — 0) — fcfl are given for several 



values as green lines. The quantity A is a measure for 
the relative deviation of the modulus of the propagation 
constant for oblique incidence from that for normal in- 
cidence. We have taken the moduli because all quanti- 
ties are complex- valued. The iso-error lines in the non- 
resonant regime (cj ;$ 150 THz) scale approximately with 
1/uj oc A. This results from the fact that the homoge- 
neous medium approximation improves with a decreas- 
ing unit cell size to wavelength ratio a/ A. For frequencies 
larger than 160 THz the structure becomes resonant and 
an abrupt change of the iso-error lines is clearly observ- 
able. 




FIG. 4: (color online) (a) Real part of the propagation con- 
stant as a function of the frequency and the angle of incidence 
where the plane of incidence is parallel to the coordinate axes, 
(b) The same as (a) with plane of incidence rotated by 7r/4, 
i.e. the TM direction. The green iso-error lines correspond to 
the relative deviations A e [0.02,0.1,0.5,2,5]% (from left to 
right) . 

Obviously in the resonant regime the MM can be only 
considered isotropic if the ratio a/ A is much less than in 
the non-resonant regime. In general the situation is iden- 
tical for both planes of incidence (TX- and TM-direction) . 
However, for larger frequencies the deviation of the cal- 
culated iso-frequency lines to the ideal ones is slightly 
smaller in TM- than in TX-direction but not significantly. 

Note that the introduction of the quantity A is only 
one option to quantify the deviation of the pertinent 
isofrequency surface from the ideal one. Other measures 
are also possible, where, for example, the deviation of the 
length of the actual wavevector compared to the length 
of the wavevector in an isotropic medium yields qualita- 
tively the same results. 

Before we proceed to discuss the results, some remarks 
concerning the investigated parameter space are in order. 
Here, the dispersion relation is not calculated up to the 
edge of the Brillouin zone, but for fixed tangential wave 
vectors which translate for propagating waves into cer- 
tain angles of incidence. This is more useful since for 
very low frequencies it is not possible to provide incident 
plane waves with tangential wave vector components k t 
that are in the order of tt/cl where a is the lattice con- 
stant. Furthermore, plane waves with wave vectors in 
the dimension of the lattice vector will always sense the 
details of the periodicity. In this case the issue of an 
optically isotropic structure is pointless. Of course, also 
for small frequencies the iso-surface are not rotationally 
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invariant but these domains are in general not accessible 
to the experiment and are negligible. Assuming the outer 
medium to have a permittivity e then the line given by 
kt/ko = \fe corresponds to grazing incidence and to the 
angular domain accessible from free space. 

Up to here we have only discussed the properties of 
the fundamental Bloch mode with the smallest imagi- 
nary part since this mode essentially dictates light prop- 
agation. While considering the coupling of an external 
field to this Bloch mode we had to suppose that only this 
mode is excited, thus restricting ourselves to a certain 
polarization state of the incident light. Nevertheless, in 
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tive index is isotropic. So as soon as the eigenvalues for 
both Bloch modes deviate from each other, the medium 
cannot be described as being effectively isotropic. This 
deviation can also serve as a measure for a meaningful 
MM homogenization. 

In Fig. [5] (a) we show the real and the imaginary part 
of the propagation constant for the two Bloch modes with 
the smallest losses at a certain angle of incidence and in 
Fig. [5] (b) their difference as function of the frequency. 
Clearly for low frequencies the propagation constants for 
both modes are identical. As soon as the frequency in- 
creases the propagation constants start to deviate from 
each other, hence the SCiC cannot be considered as op- 
tically isotropic anymore. The deviation is of course the 
stronger the larger the angle of incidence is. At normal 
incidence they are identical as explained in detail before. 
The conclusions to be drawn from these figures are of 
course in line with those from the angular dependency 
of the first Bloch mode only. The structure can be con- 
sidered as being isotropic only for small frequencies and 
a limited range of angles of incidence where this angular 
range is the larger the smaller the frequency is. 



FIG. 5: (color online) (a) Real part (solid line) and imaginary 
part (dotted line) of the propagation constant /c[/im _1 ] of the 
two lowest damped Bloch modes at k t = ko that are degen- 
erated at normal incidence, (b) difference between the real 
parts (blue line) and the imaginary parts (green line) of these 
two modes as function of the frequency. In the quasi-static 
limit, i.e. A — »• 00 the difference is of course vanishing and 
both modes are degenerated. 

the general case of arbitrarily polarized incident light one 
has to discuss both eigenstates of the polarization sepa- 
rately, as both may exhibit different propagation con- 
stants. Assuming an isotropic homogeneous medium the 
eigensolutions are always degenerated, i.e. they have the 
same propagation constant and it is possible to define 
the eigenmodes as orthogonally and linearly polarized 
plane waves. This is also the case for normally incident 
light on a C4 symmetric structure as the isotropically 
designed SCiC metamaterial. Here the light encounters 
the same physical structure at normal incidence for any 
polarization, hence the eigensolutions are degenerated. 
For oblique incidence the situation changes dramatically 
because both linear polarized eigenstates may encounter 
different structures. Both eigensolutions may have differ- 
ent eigenvalues for the propagation constant. One may 
compare the situation with an uniaxial crystal, where 
the optical axis is aligned with the z-axis. For normally 
incident light the structure can be considered as being 
isotropic, whereas at oblique incidence the eigenvalues 
for the two eigenpolarizations are different. Note that 
the eigenpolarizations are always orthogonally, linearly 
polarized for planes of incidence being identical with mir- 
ror planes of the metamaterial. 

Considering the SCiC as an effectively isotropic homo- 
geneous medium, we expect that both the two fundamen- 
tal Bloch modes are degenerated and the effective refrac- 



V. CONCLUSION 

In this contribution we have investigated the optical re- 
sponse of a highly symmetric SCiC MM in the negative 
refraction regime. The SCiC exhibits cubic symmetry, 
i.e. the highest possible symmetry for periodic meta- 
materials, and is therefore considered a very promising 
candidate for an isotropic negative refractive index ma- 
terial. By investigating at first a finite number of func- 
tional layers and the convergence of the corresponding 
effective parameters to the values obtained from the dis- 
persion relation of the infinite structure we can conclude 
that the SCiC fulfills all requirements of a homogenizable 
metamaterial. Already a single functional layer can be 
described by its bulk properties as only a single Bloch 
mode determines light propagation inside the structure. 
Nevertheless our investigation of the iso-frequency sur- 
face of the dispersion relation clearly shows that even 
this metamaterial with the highest symmetry is far away 
from having an optically isotropic response in the region 
of negative refraction. Hence, the key result of this con- 
tribution is that one must not conclude from high sym- 
metry on an optically isotropic response. Hence, optical 
isotropy of a MM can be only deduced from a careful 
inspection of the dispersion relation by taking into ac- 
count the material dispersion of metal too. This is an 
important message to the designers of isotropic optical 
metamaterials. 
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